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Abstract 

We show that two evanescently coupled x^"^^ parametric downconverters inside a Fabry-Perot 
cavity provide a tunable source of quadrature squeezed light, Einstein-Podolsky-Rosen correlations 
and quantum entanglement. Analysing the operation in the below threshold regime, we show how 
these properties can be controlled by adjusting the coupling strengths and the cavity detunings. As 
this can be implemented with integrated optics, it provides a possible route to rugged and stable 
EPR sources. 

PACS numbers: 42.50.Dv,42.65.Lm,03.65.Ud 
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I. INTRODUCTION 



The Ei„s.e.„, PodoMy and Rosen (EPR) paradox ste.ns fton. a fa^ons pape. pnbUshed 

in 1935 [1], which showed that local realism is not consistent with quantum mechanical 
completeness. A direct and feasible demonstration of the EPR paradox with continuous 
variables was first suggested using nondegenerate parametric amplification (also known as 
the OPA) P|. The optical quadrature phase amplitudes used in these proposals have the 
same mathematical properties as the position and momentum originally used by EPR. Even 
though the correlations between these are not perfect, they are still entangled sufficiently to 
allow for an inferred violation of the uncertainty principle, which is equivalent to the EPR 
pa^ado. BQ. An e.pe.n,enta, den,„n..a«cn oUM. p.„p„.a, . soon followed, 

showing a clear agreement with quantum theory ^. 

In this work, rather than using the nondegenerate optical parametric oscillator (OPO), 
we consider an alternative device using two degenerate type I downconverters inside the 
same optical cavity, and coupled by evanescent overlaps of the intracavity modes within the 
nonlinear medium. Generally, such a device may be considered as either a single nonlinear 
crystal pumped by two spatially separated lasers, or two waveguides with a x*^^^ component. 
We calculate phase-dependent correlations between the two low frequency outputs of the 
cavity in the below threshold regime, showing that this system exhibits a wide range of 
behaviour and is potentially an easily tunable source of single-mode squeezing, entangled 
states and states which exhibit the EPR paradox. The spatial separation of the output modes 
means that they do not have to be separated by optical devices before measurements can 
be made, along with the unavoidable losses which would result from this procedure. The 
entangled beams produced can be degenerate in both frequency and polarisation, unlike 
those of the nondegenerate OPO, and would exit the cavity at spatially separated locations. 
These correlations are tunable by controlling some of the operational degrees of freedom of 
the system, including the evanescent couplings between the two waveguides, the input power 
and the cavity detunings. 

The term nonlinear coupler was given to a system of two coupled waveguides without 
an optical cavity by Pefina et al. Q]. Generically, the device consists of two parallel optical 
waveguides which are coupled by an evanescent overlap of the guided modes. The quan- 
tum statistical properties of this device when the nonlinearity is of the x*^^"* type have been 
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theoretically investigated, predicting energy transfer between the waveguides Ql and the gen- 
eration of correlated squeezing j^. Coupled x^"^^ downconversion processes in the travelling 
wave configuration have also been examined theoretically, predicting that light produced in 
one of the media can be controlled by light entering the other j^, and that such a device 
can produce entanglement of the output beams The coupler with x^^^ nonlinearity 



held inside a pumped Fabry-Perot cavity, and opera ting in the second harmonic generation 
(SHG) configuration, was introduced by Bache et al. who named it the quantum optical 
dimer by analogy with various systems that display coupling between discrete sites. They 
analysed intensity correlations between the modes, predicting noise suppression in both the 
sum and the difference. 

As the intracavity x*^^^ downconversion processes have long been appreciated as sources 
of quantum states of the electromagnetic field (See Martinelli et al. for an overview), 
we will combine and extend these previous analyses to consider two coupled downconverters 
operating inside a Fabry-Perot cavity. The advantage of this proposal is the all-integrated 
nature of the device, which promises greatly increased robustness. Additional potential ad- 
vantages are the reductions in threshold pump power and phase noise, relative to current 
practise. Another potential advantage as compared to the normal type II polarisation nonde- 
generate OPO lies in the difficulty of stabilising this device at frequency degeneracy 
Our proposal should be well stabilised by the linear coupling, without having to add any 
additional features. 



II. THE SYSTEM AND EQUATIONS OF MOTION 

The physical device we wish to examine differs from that described in Ref. ^| in one 
important detail. We will analyse it in the downconversion regime, where the cavitv pumping 
is at a frequency 2^^ ~ cj^,. As this device has been described in detail in Ref. jllj, we will 
give a briefer description of the essential features here. The system of interest consists 
of two coupled nonlinear x*^^'' waveguides inside a driven optical cavity, which may utilise 
integrated Bragg reflection for compactness. Each waveguide supports two resonant modes at 
frequencies oJa, i^b, where 2uJa — oJh- The higher frequency modes at ujb are driven coherently 
with an external laser, while the nonlinear interaction within the waveguides produces pairs 
of downconverted photons with frequency uOa ■ We assume that only the cavity modes at 
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these two frequencies are important and that there is perfect phase matching inside the 
media. The two waveguides are evanescently coupled as in Ref. [uj. Besides the differences 
in the pumping frequency, we will be interested in the phase-dependent correlations necessary 
to demonstrate entanglement and the EPR paradox, rather than the intensity correlations 
of Ref. 

The effective Hamiltonian for the system can be written as 

"^e// 'Hint ~l~ 'Hcouple ~l~ 'Hpump ~l~ T~(-resi (1) 

where the nonlinear interaction with the x*^^"* media is described by 

a\ %i - djbl + al %2 - ajbl] . (2) 



Hint = ifi— 



Here k denotes the effective nonlinearity of the waveguides (we assume that the two are 
equal), and a^, hk are the bosonic annihilation operators for quanta at the frequencies 
uja, UJ), within the crystal k (= 1,2). The coupling by evanescent waves is described by 



H-couple 



aia\ + a[a2 



+ hJb 



bibl + blk] , (3) 



where the Jk are the coupling parameters at the two frequencies, as described in Ref. 
We note that in that work it is stated that the lower frequency coupling, J^, is generally 
stronger than the higher frequency coupling, Jf,, and also that values of Ja as high as 50 
times the lower frequency cavity loss rate were calculated to be physically reasonable. The 
cavity pumping is described by 



tib\ - e*6i + e2b\ - e\b^ , (4) 



where the represent pump fields which we will describe classically. Finally, the cavity 
damping is described by 

2 

Ke. = ^$^(r^4 + r'^l)+^-c-, (5) 

where the represent bath operators at the two frequencies and we have made the usual 
zero temperature approximation for the reservoirs. 

With the standard methods jl^, and using the operator/c-number correspondences (hj ^ 
a, , 6, ^ /3i ), the Hamiltonian can be mapped onto a Fokker-Planck equation for the Glauber- 

n n 

Sudarshan P-distribution [ly, li2|. However, as the diffusion matrix of this Fokker-Planck 



equation is not positive-definite, it cannot be mapped onto a set of stochastic differential 
equations. Hence we will use the positive-P representation [3| which, by doubling the 
dimensionality of the phase-space, allows a Fokker-Planck equation with a positive-definite 
diffusion matrix to be found and thus a mapping onto stochastic differential equations. 
Making the correspondence between the set of operators {aj, a], bj, Sj) (j = 1, 2) and the set 
of c-number variables {aj,aj' , Pj, we find the following set of equations, 

^ = -(7a + + Ka|/?i +iJaa2 + a/k^'7i(^)> 

da^ 

HT 
da2 
~dt 
da2 

~dr 

— = ei - (7i, + %l\h)p\ - -jOix + %JhP2-, 

^ = - (7. - .A,)/5+ ^Ji^^t, 
dp2 

dt ' 2 

dt 



-(7a + iAa)a2 + Ka^p2 + iJa^l + \f^2 Vsit): 



-(7a - iAa)a^ + Ka2P2 - '^■Ja^t + 



£2 - {lb + ^Ab)/?2 - + iJbl3i, 

el - (7fe - i^b)Pt - S «^ ' - ^JbPt. (6) 



where the 7^ represent cavity damping. We have also added cavity detunings Aa,?, from the 
two resonances, so that for a pump laser at angular frequency 2ujlj one has Aa = uJa — oJiajnd 
A5 = ojh—^ojL- Below, in sectionEl we will investigate detuning effects in greater detail. The 



real Gaussian noise terms have the correlations rij{t) = and rij{t)rik{t') = 6jkS{t — t'). Note 
that, due to the independence of the noise sources, (Pk) and (/5^) are not complex 
conjugate pairs, except in the mean over a large number of stochastic integrations of the 
above equations. However, these equations do allow us to calculate the expectation values 
of any desired time-normally ordered operator moments, exactly as required to calculate 
spectral correlations. 

III. LINEARISED ANALYSIS 

In an operating regions where it is valid, a linearised fiuctuation analysis provides a simple 
way of calculating both intracavity and output spectra of the system (l^,l20|, by treating it as 
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an Ornstein-Uhlenbeck process |21|. To perform this analysis we first divide the variables of 
Eq.Elinto a steady-state mean value and a fluctuation part, e.g. ai af' + Sai and so on for 
the other variables. We find the steady state solutions by solving the equations ^ without 
the noise terms (note that in this section we will treat all fields as being at resonance), and 
write the equations for the fiuctuation vector 6x = [6ai, 6af, 6a2, Sa2, SPi, 6P2, ^PtV^ 
to first order in these fiuctuations, as 



d6x = -A6x dt + BdW, 



(7) 



where the drift matrix is 



with 



A., 



A, 



A, 



A 



Aaa 
Aba 



- A* 



A 



bb 



la 




-iJa 





-K(3i 


la 





iJa 


-iJa 


la 







iJa 




la 















nal'* 


















Kan 















lb 

76 
-iJb 

iJb 



-iJb 

iJb 

lb 

lb 



(8) 



(9) 



In this equation, dW is a vector of real Wiener increments, and the matrix B is zero except 
for the first four diagonal elements, which are respectively , ^/k^, ^/k^^. 

The essential conditions for this expansion to be valid are that moments of the fluctuations 
be smaller than the equivalent moments of the mean values, and that the fluctuations stay 
small. In the case of the single optical parametric oscillator (OPO), it is well known that 
there is a critical operating point around which this condition does not hold. This point is 
easily found by examination of the eigenvalues of the equivalent fluctuation drift matrix for 
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that system, and this procedure is also valid in the present case. The fluctuations will not 
tend to grow as long as none of the eigenvalues of the matrix A develop a negative real part. 
At the point at which this happens the linearised fluctuation analysis is no longer valid, as 
the fluctuations can then grow exponentially and the necessary conditions for linearisation 
are no longer fulfilled. In this work we will only be interested in a region where linearisation 
is valid. 

To examine the stability of the system, we first need to find the steady state solutions for 
the amplitudes, by solving for the steady state of Eq. El with the noise terms dropped. As in 
the usual optical parametric oscillator, there is an oscillation threshold below which a^"* = 
and only the high frequency mode is populated. In the present case, for a real pump, we 
find Z?!** = e/(7f, — iJb)- Inserting these solutions in the matrix A allows us to find simple 
expressions for the eigenvalues, 

Ai,2 = Ib + iJb, 
^3,4 = lb- iJb, 

A7,8 = 7a - ^J[>^'ey^!-J!]. (10) 

Here we have introduced auxiliary variables, ^a,b = \Jl1b~^ ^ab- immediately see that 
Ay^s can develop negative real parts for a pump amplitude greater than the critical value, 
Cc = laHb/ K'- As it must, this expression reduces to the single OPO expression of •jalb/i^ when 
the couplings are set to zero. In that case, there is then a stable above threshold solution 
in which the high frequency mode inside the cavity remains constant, independently of any 
further increase in the pumping, and the low frequency mode becomes occupied. 

In the present case, it is not simple to find general expressions for these above threshold 
solutions analytically, but as we will concentrate our attention on the rich variety of below 
threshold behaviour which is exhibited, this is not important here. We note here that, unlike 
the single OPO case with a resonant cavity, the threshold pumping is not a constant for fixed 
cavity loss rates, but is a function of the coupling strengths between the two waveguides. 
Using the below threshold solutions, we may calculate any desired time normally-ordered 
spectral correlations inside the cavity using the simple formula 

S{uj) = {A + tool)-^ BB"^ {A^ - tooiy^ , (11) 
7 



after which we use the standard input-output relations to relate these to quantities 
which may be measured outside the cavity. 



IV. QUANTUM CORRELATIONS 
A. Single mode squeezing 

The first quantities we wish to calculate are the single mode quadrature squeezing spectra, 
to compare these with the well-known results for the normal uncoupled OPO. Defining the 
quadrature amplitudes as 

X| = ajQ-'' + a]e'\ (12) 
(where j = 1, 2), we will use the notation 

X° = Xj, 

XJ = F,. (13) 

We note here that the quadrature definitions do not need to specify whether it is mode a 
or b which is involved, as we do not find any interesting behaviour in the high frequency 
modes below threshold and hence will only present results for the low frequency modes. 
With this normalisation the coherent state value for the quadrature variances is one. To 
simplify our results we will assume that the pumping terms for each crystal are real and 
equal (ei = e2 = e). 



The expressions for the be 



OPO case are well known 



I22I 



ow threshold low frequency quadrature variances in the single 
being 



',.A - 1 ^ 



X — J- -r 7 r^— 29' 

[lalb - Ke)2 + 7^^^2 

7» = I . ..... ^"°r^.... . (14) 



and predicting zero-frequency squeezing which becomes perfect in the Y quadrature as the 
pump approaches the critical threshold value, e = '~ialb/K', although the linearised analysis 
breaks down near this point. Note that the variances inside and outside the cavity are 
related by 5^"* = l + 27aF(X). Our coupled system would be expected to exhibit the above 
values in the limit as Ja^b 0, which provides a standard for comparison with the analytical 
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results. In the general case, we find that S'^g = 5*°"*, as expected. We also find that the 
coupling means that the intracavity high frequency field is no longer real, but has a phase 
given by 6fe = tan"^( Jfe/76). 

This will mean that the optimum correlations will no longer generally be found in the Xj 
and Yj quadratures, but at some other phase angle, as found previously for second harmonic 
generation in detuned cavities 0|. Experimentally, this does not present a problem as the 
local oscillator phase is normally swept across all angles, which must therefore include the 
optimum angle. We can find analytical solutions for the angle of maximum single-mode 
squeezing (and antisqueezing) , for example, these differing by 7r/2 and being found as 



Oo^t = tan-i { ^ 2V{X,Y) ^ ^^^^ 



V{Y)-ViX)±^^ViY)-ViX)^ +4[V(X,F) 

where V{A,B) = {AB) — {A){B). However, as this expression is a complicated function of 
several variables when written out in full, and will not necessarily give the optimum choices 
at all frequencies, nor when we consider correlations between the modes, we will present 
results where the local oscillator angle has been optimised numerically. 
The X and Y spectral variances outside the cavity are found as 

^out ( ) = 1 , 47,^e {7fc [7fc^(u;^ - J^J + J^l + (7a7b + ^^f\ + '^laJl^^) 
^outf ^ _ 1 47.^6 {7b [%\uj' - J^) + J^-fl + (7a7fc - i^ef] + '^'jaJ^^e] 

47^6 ^ + [76 (7a - - K2e2] 

which, as expected, reduce to the single OPO expressions above (fT^ when the coupling 
terms are set to zero. The output covariance is 

^ ~ .2-4 2^ r~2.~2 7 ^ ^ 

47a7b ^ + [lb [la - UJ^) - K^e^\ 

which will give Oopt = 0,7r/2 for the uncoupled case, where Y is the squeezed quadrature 
and X the antisqueezed quadrature. 

In Fig. [U we show the single-mode output spectral quadrature variances for the quadrature 
of best squeezing as the low-frequency mode coupling strength is varied, beginning with 
Ja = Jb = la = lb = 1- Wc uotc hcrc that the pump values used in all the displayed results, 
ej = 0.5ec, depend on the couplings as stated above and are therefore different for different 




Figure 1: S'^g{io) for 7 = 1, = 1 and different Ja, all at the 9 of maximum squeezing. The 
solid line is for = 1 and 9 = 113°, the dash-dotted line is Ja = 2, the dotted line is = 5 and 
the dashed Une is Ja = 10, all for 9 = 22". The pump amplitude is e = 0.5ec in each case and all 
quantities plotted in this and subsequent graphics are dimensionless. Note that all plotted spectra 
are symmetric about zero frequency and all results shown use the value k = 0.01 and 7a = 76 = 7- 

combinations of the couplings, but are all the same fraction of the threshold value. We find 
less single-mode squeezing than in the uncoupled case for the same ratio e/ec, and also find 
that changing Jh mainly serves to change the angle of maximum squeezing. Changing Ja 
changes the frequency at which the maximum of squeezing is found. We see that this device 
is not as efficient at producing squeezed single-mode outputs as the normal OPO, but as 
we are interested in the quantum correlations between the two output modes, we will now 
examine these. 
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B. Entanglement and the EPR paradox 



An entanglement criterion for optical quadratures has been outlined by Dechoum et 
al. following from criteria developed by Duan et al. [3| which are based on the insepa- 



rability of the system density matrix. A theoretical met 



lod to demonstrate the EPR paradox 



using quadrature amplitudes was developed by Reid ^, using the mathematical similari- 
ties of the quadrature operators to the original position and momentum operators. We will 
briefly outline these criteria here and then apply them to our system, using the quadrature 
operators Xj and Yj. Note that even though these quadratures have the same mathematical 
properties as the canonical position and momentum operators for the harmonic oscillator, 
they correspond physically to the real and imaginary parts of the electromagnetic field, not 
its position and momentum. 

To demonstrate entanglement between the modes, we define the combined quadratures 
X± = Xi ± X2 and Y± = Yi ± Y2 and calculate the variances in these, which we may do 
analytically. Optimising the resu 
Following the treatment of Ref. 



t for arbitrary phase angles is better performed numerically. 
. 3|, entanglement is guaranteed provided that 



S°j^l + S^^ < 4. (18) 

We note here that the combined variance defined in this way has an obvious relationship 
with the well-known two-mode squeezed states which are produced, for example, by the 
nondegenerate OPO 0, 0|i but that the quadratures between which we find entanglement 
here are not the same as those which are entangled in that case, where these are X_ and Y^. 
In the present case, considering only the phase angles 9 = and tt/2, we find entanglement 
with and YL. The two individual variances can be written as 

gout ^ gout ^ gout ^ 2T/(Xi, X2), 

Sxl = ^x' + S"x^^ ± 2\/(fi, (19) 

The individual quadrature variances are given above (fT6|) . while for the covariances we find: 

V{X„ X,) = -^JMll-e _ ^20) 

and V{Yi, Y2) = —V{Xi, X2), showing that the X quadratures are anticorrelated and the Y 
quadratures are correlated. Although these results allow us to write analytical expressions 

11 



for the combined variances, these are rather bulky and not very enlightening, so we will not 
reproduce them here. 

To optimise the degree of entanglement as a function of the quadrature phase angle, we 
investigate the output spectral correlation 

sr{x-)+sr{Y+), (21) 

where the X quadratures are at the angle 9 and the Y quadratures at the angle 6 + 71/2. 
What we find, as shown in Fig. [2l is that the degree of entanglement and the frequency at 
which it exists depend on the coupling strength Ja while the optimum angle depends on Jb. 
When we hold Ja constant and increase Jb, we find that the maximum of entanglement is 
always found at zero frequency, but that the optimum quadrature angle changes. 

To examine the utility of the system for the production of states which exhibit the EPR 



paradox, we follow the approach of Reid [Sj- We assume that a measurement of the Xi 
quadrature, for example, will allow us to infer, with some error, the value of the X2 quadra- 
ture, and similarly for the Yj quadratures. This allows us to make linear estimates of the 
quadrature variances, which are then minimised to give the inferred output variances, 

^ 1 2 



(Xi) = 5°"* 



V{X^,X2) 



X2 

V{Yi,Y2)Y 

S^^^Yi) = - ^ oout ^ • (22) 

The inferred variances for the j = 2 quadratures are simply found by swapping the indices 
1 and 2. As the Xj and Yj operators do not commute, the products of the variances obey 
a Heisenberg uncertainty relation, with S'^^S^* > 1. Hence we find a demonstration of the 
EPR paradox whenever 

-5°^ (x,)5°f (y,) < 1. (23) 

With the expressions for the variances given in Eq.[Tniand the covariances of Eq.lSHl we have 
all we need to calculate the EPR correlations. Once again, however, the full expressions are 
somewhat unwieldy, so we will present the results graphically. 

In Fig. ini we present the results for optimised quadrature phase angles while Jb is held 
constant at a value of 7 while Ja is increased. Note that again the angle 9 refers to the X^ 
quadratures, while the conjugate quadratures are at an angle of ^ + 7r/2. Changing Jb serves 
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1 1 ^ ^ ^ 1 

5 10 15 20 

CO (units of y) 

Figure 2: Demonstration of entanglement, using S^^^{X^) + S^^^{Y+), for 7 = 1, = 1, and Ja = 1 
(solid line), 2 (dash-dotted line), 5 (dotted line) and 10 (dashed line). The quadrature angle for X 
is 67° and that for Y is 157°. The pump amplitude is e = O.Sec- 

to change the angle of the maximum violation, without changing the degree of violation, 
while changing Ja changes both the degree and the frequency of the maximum violation. As 
expected, these results are the same for both outputs of the device. 

V. DETUNING THE CAVITY 

Often in optical systems the best performance is found when the cavity is resonant for the 
different modes involved in the interactions. In the present case we find that detuning the 
cavity by the appropriate amount from the two frequencies allows for some simplification 
of the theoretical analysis and can actually improve some quantum correlations. With 
detunings included, the steady state below threshold solutions for the high frequency mode 
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1.4 




Figure 3: Demonstration of the EPR correlation for Jf, = 7 = 1 and Ja = 1 (solid line, 9 = 67°), 2 
(dash-dotted line, d = 67°), 5 (dotted line, 6 = 113°) and 10 (dashed line, 9 = 113°). The pump 
amplitude is e = 0.5ec. 

are found as 

[7b - H-'6 - ^b)\ 

so that, setting A;, = J;,, we return to the well-known real solutions for a single OPO. If we 
then set = J^, define the new variables Ap = ai + a2 and = ai — 02, and eliminate 
the time dependence of Pa,b, we can write positive-P stochastic equations as 

^ = -'JaAp + KpssA'^ + X^KpTs {ill + ^3) , 

d-A^ I 

= -7a^p + I^PssAp + V I^Pss (,V2 + Vi) 5 

= - [7a + 2i Ja] A^ + KpssA'^ + {iji - 7]^) , 

= - [-fa- 2iJa] A+ + KpssAm + {V2 ' V^) ■ (25) 
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In the above, the noise terms are the same as those of Eq. El We note here that, although 
it is the detuning in the low frequency mode that allows us to write the equations for Ap 
and Ap in a particularly simple form, Af, also plays a role in that it allows us to treat Pss 
as real, which will make the interesting quantum correlations in and between the X and Y 
quadratures, so that we do not have to examine all possible local oscillator angles to find 
the best performance. 

Following the same linearisation procedure as in section llllj we find the corresponding 
drift and noise matrices. 



A. 



pm 
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pm 
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In terms of the quadratures used in section IIV[ we now define 



(26) 



(27) 



Xp — Ap + Ap — Xi + X2, 
Xm = Am + A'^ = Xi — X2, 
Yp = -I [Ap - A+) = Fi + F2, 
Ym = -% {Am - A+) = Fi - 1^2, 



(28) 



and give expressions for the output spectral variances of these new quadratures. For the Xp 
and Yp quadratures these are particularly simple, 



sT[^) = 2 + 



ilalb-^e) +lb^'^' 

8lalbfie 



(29) 



ilalb + ne) + iluj"^ 

and are readily seen to be the sum of the variances for two uncoupled OPOs, as given in 
Eq. [m As in that case, the zero-frequency variance in Yp is predicted to vanish at the 
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critical pump value of ec = "^alb/ 1^, although, as should be well known, a linearised analysis 
is not valid in that region. However, the degree of squeezing is more than was found to 
be available in the doubly resonant case considered above. The other two variances do not 
uncouple and have more complicated expressions, 

Qout/ X ^ o , [(7,76 + «;e)2- 7^(472 -^2)] 



STiuj) 



_ r, 87,7b^e [(7,7fe - Kef - 7,^(47,^ - u^)] 
'"^^ + AJl - co^) - K^e^f + A^l^tu^ ■ 

(30) 

Graphical results for the combined quadratures which exhibit squeezing are shown in 
Fig. m from which it is obvious that by far the best squeezing quadrature is 1^, which, for 
these parameters, shows almost 90% squeezing at zero-frequency. The quadratures and 
Ym show only a very small degree of squeezing far from zero frequency. What this result 
shows, along with the results for the resonant cavity, is that the low frequency modes want to 
oscillate at two distinct frequencies, as is normal for coupled systems. The detuning chosen, 
Aa = J a, moves the sum mode frequency closer to resonance while the other frequency is 
further detuned. Along with the choice of A;, so as to make the intracavity high frequency 
amplitude real, this results in maximised single-mode noise supression and entanglement 
centred on zero frequency. 

Using these results, we can now investigate the degree of entanglement, as done above 
for the resonant case. As shown in Fig. El we find that the quadratures Yp and Xm are 
entangled, exactly as in the single OPO case. As with the squeezing, the detunings have 
moved the maximum of entanglement to zero frequency. A sign of the out of resonance 
mode attempting to resonate is seen in the small degree of entanglement apparent around 
uj ^ 2O7. We also see that the degree of entanglement is less than in the case with zero 
detuning, shown previously in Fig. O although it must be remembered that the absolute 
pump powers are not the same, merely the ratios e/ec- Finding analytical expressions for 
EPR correlations is not possible using this coupled-mode approach, as, although we can 
calculate the necessary covariances, for example, V"(Xi,X2) = [V^(Xp) — V^(Xm))]/4, it is 
not obvious how to separate out the single-mode variances. However, these can still be 
calculated numerically using the full single-mode equations with the appropriate detunings. 
That the system clearly demonstrates the EPR paradox is shown in Fig. El although again 
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Ql > > > > 1 

5 10 15 20 25 

(0 (units of y) 

Figure 4: Output spectral variances of the combined modes for J;, = A;, = 7^ = 7^ = 7 = 1 and 
J„ = Aa = 10. The sohd hne is S^f, the dash-dotted hue is and the dotted hue is S^^. A 

value of less than 2 represents squeezing. The pump amplitude is e = 0.5ec. 

the maximum inferred violation is less than in the resonant case. 

We note here that all the quantities shown for the detuned system are actually calculated 
at a lower absolute pump power than in the resonant case. For positive detunings, the 
critical pump amplitude is found as 

ec = ^bl + iJa-Aam! + U-A,y]/K , (31) 

SO that our choice of detunings means this is no longer a function of the coupling strengths. 
Therefore a careful choice of detunings has two main advantages in that it fixes the quadra- 
tures for which the maxima of quantum features are found, and means that the pumping 
necessary to a good performance does not vary with the coupling strengths. 

The choice of detunings shown has the possible disadvantage that, as the effective coupling 
is now only in the mode, which is moved away from resonance, the quantum correlations 
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Figure 5: The sum of the output spectral variances, + Sy"*, for = A;, = 7a = 7b = 7 = 1 
and Ja = Aa = 1 (solid line), 10 (dash-dotted line) and 20 (dashed line). A value of less than 4 
represents entanglement. The pump amplitude is e = 0.5ec. 

which depend on both the modes are slightly degraded. This is readily seen from the figures 
because those correlations which include and Y^n change more with Ja than do the others. 

VI. CONCLUSION 

This system exhibits a wide range of behaviour and is potentially an easily tunable source 
of single-mode squeezing, entangled states and states which exhibit the EPR paradox. The 
spatial separation of the output modes means that they do not have to be separated by 
optical devices before measurements can be made, along with the unavoidable losses which 
would result from this procedure. The entangled beams produced can be degenerate in 
both frequency and polarisation, unlike those of the nondegenerate OPO, and would exit 
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Figure 6: The product of the inferred output spectral variances, S^^'^ Sy^ , for Jf, = Af, = 7^ = 7b = 
7 = 1 and Ja = Aq = 1 (solid line) and 10 (dash-dotted line). On this scale, the result for Ja = 20 
is indistinguishable from that for Ja = 10. A value of less than 1 represents a demonstration of the 
EPR paradox. The pump amplitude is e = O.Scc. 



the cavity at spatially separated locations. This may be a real operational advantage over 
the nondegenerate OPO, which is also known to produce nonclassical states. The tunability 
that exists because of the number of different parameters which can be experimentally 
accessed, such as the coupling strength, the pump intensity and the detunings, may make 
it interesting for a range of potential applications which would require the availability of 
states of the electromagnetic field with varying degrees of nonclassicality. Since this type of 
system is compatible with integrated optics techniques, it may provide a more robust source 
of entanglement than interferometers that use discrete optical components. 
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